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Abstract 

We show that for a general four derivative theory of gravity, only the holographic 
entanglement entropy functionals obey the second law at linearized order in perturbations. 
We also derive bounds on the higher curvature couplings in several examples, demanding 
the validity of the second law for higher order perturbations. For the five dimensional 
Gauss-Bonnet theory in the context of AdS/CFT, the bound arising from black branes 
coincides with there being no sound channel instability close to the horizon. Repeating the 
analysis for topological black holes, the bound coincides with the tensor channel causality 
constraint (which is responsible for the viscosity bound). Furthermore, we show how to 
recover the holographic c-theorems in higher curvature theories from similar considerations 
based on the Raychaudhuri equation. 
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1 Introduction 

In the context of AdS/CFT, entanglement entropy of the boundary quantum field theory can 
be calculated using the Ryu-Takayanagi DEI prescription and its generalizations to higher 
curvature gravity theories. The derivation of this holographic entanglement entropy for two 
derivative Einstein gravity was proposed by Lewkowycz and Maldacena in [3j. This derivation 
allows one to obtain a surface equation (for static situations) for the entangling surface which 
minimizes an entropy functional. There were attempts to extend this calculation to find the 
entangling surface equation in higher curvature theories in [4] and the corresponding entropy 
functionals were argued to be different from the Wald entropy PS S3 d- In particular, for 
Lovelock theories the entropy functional is the so-called Jacobson-Myers (JM) [8, 9] one while 
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for general four derivative theories, it coincides with the Fursaev-Patrushev-Solodukhin (FPS) 
entropy functional ra. However, there still exist problems in finding the entangling surface 
using the Lewkowycz-Maldacena method and there are potential ambiguities related to higher 
order extrinsic curvature terms in the entropy functionals. Moreover, it leads one to wonder if 
such entropy functionals could arise independent of AdS/CFT from different and perhaps more 
fundamental considerations. 

Wald and collaborators mm had established the first law for black hole mechanics for 
any diffeomorphism invariant theory of gravity and proposed that the entropy of a black hole 
is a Noether charge associated with the Killing isometry generating the horizon. The Wald 
entropy suffers from various ambiguities [13] and therefore does not provide a unique answer 
for the horizon entropy. Although, none of these ambiguities contribute to black hole entropy 
of a stationary Killing horizon with regular bifurcation surface, the study of the second law 
of black hole mechanics beyond general relativity (GR) shows that these ambiguities need to 
be carefully included in the expression of black hole entropy to obtain an increase theorem 
similar to Hawking area theorem in GR. In fact, for black holes in Lovelock gravity, it is the 
JM functional which leads to an increase theorem for linearized perturbations [IT] , 

Recently, it has been pointed out that the 2nd law for spherically symmetric black holes 
is satisfied at linear order in perturbations in general four derivative theories of gravity by 
the holographic entanglement entropy (HEE) functionals j!5| . This has been also generalized 
beyond spherical symmetry and to generic perturbations up to linear order [16] and a detailed 
construction of the entropy functionals has been proposed for any metric theory of gravity 
such that the linearized second law holds true. For theories with higher curvature terms, the 
construction produces the HEE functionals. This is a remarkable result since it allows for a 
way independent of AdS/CFT of deriving entropy functionals-these same entropy functionals 
therefore find applications in diverse settings. However, what still needs to be addressed is if 
only the HEE functionals alone do this job. In this paper, we first show that at linear order in 
perturbations only the HEE functionals can satisfy the second law. We start with an expression 
for horizon entropy with all possible ambiguous terms relevant for linear order in perturbation. 
We fix these terms by demanding the validity of the linearized second law and the final entropy 
coincides with HEE functionals. To be clear, since our analysis is going to be perturbative, 
we will not be able to rule out higher order extrinsic curvature terms in the functionals (e.g. 
0(/F 4 ) in the FPS functional)-see section 5. 

We will tackle the question keeping in mind two different motivations. First, we are in¬ 
terested in what happens in the case of asymptotically flat black holes. We will consider two 
different types of perturbations-the first kind will be due to radially symmetric, slowly falling 
matter and the second kind being a shear perturbation. We will find that for the Gauss-Bonnet 
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(GB) theory, if we go to second order in perturbations, the second law is automatically satis¬ 
fied for regular horizons if the GB coupling is positive. If the GB coupling is negative then a 
certain lower bound on the horizon radius ensures that the second law is satisfied; since this 
suggests that such black holes cannot be formed from collapsing matter, one can take this as 
an indication that the negative coupling is disfavoured. If we consider a Ricci-square theory 
then we find that similar extra conditions on the mass of the black hole may be needed for the 
second law to hold. Second, we are interested in what happens in the context of AdS/CFT. 
We will find that for the GB theory, the second order analysis leads to a bound on the Gauss- 
Bonnet coupling. For black branes, this bound coincides with the absence of sound channel 
instabilities at the black hole horizon. A similar analysis for the Ricci- square theory bounds 
the corresponding coupling constant. For GB topological black holes im, for the zero mass 
case we find that the second law bound coincides with the tensor channel causality constraint 
[IS, 19]. This bound is what leads to the lower bound on the ratio of shear viscosity to entropy 
density in the dual plasma and agrees with the lower bound on the a/c ratio in 4d CFTs 1231 
with J\T = 1 supersymmetry. We also investigate quadratic theories with and R 2 terms. 

In these theories, the causality constraints following from [20] are trivial whereas the second 
law bound in certain examples we study is nontrivial. 

Using the HEE entropy functionals and the Raychaudhuri equation, we then turn to the 
question of holographic c-theorems. We will argue that the holographic c-functions found in 
Eza arise naturally from such a consideration provided the matter sector satisfies the null 
energy condition. In [22] . a derivation for the c-functions in Einstein gravity, found in [23], was 
given using the Raychaudhuri equation. In [24] an attempt was made to extend this to higher 
curvature theories using the Iyer-Wald prescription. Unfortunately the resulting functions do 
not give the correct central charge (namely the A-type anomaly coefficient in even dimensions). 
In [23J it was observed that in all sensible holographic models, which can be constructed by 
demanding the absence of higher derivative terms in the radial direction (to ameliorate the 
problem of ghosts) allow for a simple c-f nn otion which at the fixed points coincide with the 
A-type anomaly coefficient in even dimensions. This c-function was monotonic under RG flow 
provided the matter sector satisfied the null energy condition. We will find that using the HEE 
entropy functionals and considering the Raychaudhuri equation naturally leads to the same 
c-functions as in ED- 

The organization of the paper as follows: In the next section we briefly describe the am¬ 
biguities in Wald’s Noether charge construction. In the section [3] set up for proving second 
law has been described. Next, we determine the coefficients of the ambiguous terms in Wald’s 
construction using linearized second law. In section [5j we go beyond linear order in perturba¬ 
tion to study the GB theory and determine bound on the coupling parameter using AdS black 
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hole solution with different horizon topologies as background. We also consider various other 
examples like critical gravity theories and put bounds on the couplings in those theories. In 
section 6, we derive holographic c-functions using the HEE functionals and the Raychaudhuri 
equation. We end this paper by discussing several follow-up questions. We use { — 
signature and set G = 1 throughout the paper. 

2 Ambiguities in Noether charge method for non-stationary 
horizons 

Let us start by describing the geometry of the horizon of a stationary black hole in D dimensions. 

The event horizon is a null hyper-surface T~i parameterized by a non-affine parameter t. The 
vector field k a = ( d t ) a is tangent to the horizon and obeys non-affine geodesic equation V a k b = 
k k b with k is the surface gravity of the horizon (a, b, ... are bulk indices). All t = constant slices 
are co-dimension 2 space-like surfaces and foliate the horizon. We construct another auxiliary 
null normal to the t = constant slices l a , and the inner product between k a and l a satisfies 
k a l a = — 1. The induced metric on any t = constant slice of the horizon is now constructed as 
h a b = g a b + 2 k(Jb)- The horizon binomials are then given by e a b = ( k a lb — kbl a )■ 

The change in the induced metric along these two null directions can be expressed as a 
sum of trace part and a trace-free symmetric part (assuming null congruences are hypersurface 
orthogonal therefore have zero twist). The trace part measures the rate of change of the area 
of the horizon cross section along the null generators which is known to be the expansion and 
the other part measures the shear of the null geodesic congruences. We denote 9k and 9[ to be 
the two expansion parameters in these two null directions and similarly we have two shears 
and <t“ 6 . Note that 9k and cr£ 6 typically vanish on a stationary horizon but 9i and crf fe in general 
do not vanish. 

As mentioned in the introduction, the entropy of a stationary black hole in any general 
covariant theory is given by the Noether charge associated with the boost symmetry generating 
Killing vectors at the horizon mini. The Wald entropy functional for any general covariant 
Lagrangian C in D dimension takes the form: 



( 1 ) 


where C denotes any horizon slice and \fh is the area element. R a bcd is the Riemann tensor 
of bulk geometry and dA is the area of the D — 2 dimensional cross-sections of horizon. This 


formula gives a unique expression for entropy so long one is confined to any stationary slice 
or to the bifurcation surface of the horizon. However, as pointed out in |T3) [L2], the entropy 
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expression constructed via Noether charge approach has several ambiguities if one tries to apply 
it for nonstationary slices of the horizon. In fact it turns out, the Wald entropy formula is just 
one of several possible candidates for the entropy. In [13], Jacobson, Kang and Myers (JKM) 
identified three different types of ambiguities which may alter the Wald construction. Among 
these, the only relevant one for our purpose will be the one which gives us the freedom to add 
to the Wald entropy a term of the form with arbitrary coefficients: 


S 


(JKM) 


= -2tt / X.YVhdA 


( 2 ) 


where the integrand is invariant under a Lorentz boost in a plane orthogonal to an arbitrary 
horizon slice C but the terms A" and Y are not separately boost invariant [ 161 l25 j. These 
ambiguities vanish for Killing horizons but do not necessarily vanish on a nonstationary slice of 
the horizon. As a result the entropy functional for a nonstationary horizon may be expressed 

a43 


5 = 


-27r J VhdA 


dL 


dR 


-tabbed, pOkOl C[(Jk.(Jl 


abed 


(3) 


Here we have used the notation cr^erf 6 = OfcCp. As can be seen from (J 3 ]) that the ambiguous 
terms in the entropy formula involve equal number of k and l subscripts which follows from 
the fact that these are the only boost invariant combinations that can appear in the entropy 
functional. Also on a stationary slice the ambiguity terms vanish and then the expression 
coincides with the Wald formula Q. The First Law of black hole mechanics PI [13] doesn’t 
fix the coefficients p and q uniquely. So, to fix the coefficients of these terms one is thus forced 
to examine whether S obeys a local increase law. 

Recently in [16], the second law has been shown to hold for any arbitrary higher curvature 
theories if one allows only linear perturbations to a stationary black hole. This analysis also 
shows, the entropy functional proposed in the context of Holographic Entanglement Entropy 
(HEE) [5| [6, TO] remarkably matches with the S. However, it is not apparent from this anal- 
ysif0 how just one condition (entropy increases along t) fixes two coefficients in the entropy 
functional! In this paper, we will show explicitly how this fixing happens for curvature squared 
gravity theories. Note that the unknown coefficients are sitting in front of products of shear 
and expansion terms. Since shear and expansion belong to different irreducible parts of a ten- 

1 In the context of entanglement entropy, one typically evaluates this formula on a codimension-2 surface 
with one timelike nA' ,a and one spacelike rd T,a normal such that , rh 1 )® = k Jd and rk 2 )® = k ^ and the 

extrinsic curvature for that codimension-2 surface with a induced metric h a b defined as, /C^ = ^£ n (i)g a b- 
Trace of this is defined as XD = /C^h ab . Also we have = —2 dkOi and IC^ ab IC^ ab = + v < k&iai?j , 

where D is the dimension of the bulk spacetime. 

2 The criterion in [16] is sufficient but may not be necessary. 
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sor, they will remain separated even when perturbation is turned on at the linearized level. 
Consequently we will obtain two independent equations when we demand that entropy increase 
law holds for every slice of the horizon. This fixes all the coefficients of the entropy functional 
uniquely for generic quadratic curvature gravity. 

It is important to mention that we have only considered ambiguities which are quadratic 
order in expansion and shear. We may also have higher powers of such products added to 
the entropy functional. However, we cannot fix such terms using linearized second law but in 
curvature squared theories these terms do not enter at linear order ra. Our analysis matches 
with the result obtained in [T5] where the coefficients were fixed for Ricci 2 theory using a Vaidya 
like solution in the linearized increase law. 


3 Second law set up 

Let us turn to the apparatus needed to verify the second law of black hole mechanics. The 
equation of motion for a generic higher curvature metric theory of gravity will be of the form, 


Gab + H a b — 8nT a b , 


(4) 


where G a b is the Einstein tensor coming from the Einstein-Hilbert part of the action and H a b 
is the part coming from higher curvature terms in theories with a cosmological constant there 
will also be an additional term proportional to the metric which we can absorb into G a b■ T a b is 
the energy momentum tensor which we will assume to obey the Null Energy condition (NEC): 
T a bk a k b > 0 for some null vector k a . 

We will use the Raychaudhuri equation for null geodesic congruence which describes the evo¬ 
lution of the expansion along the horizon generating parameter t. In nonaffine parametrization 
this looks likcE] 

rW, Hi 

Rkk • (5) 


dd k 

—rr = K&k 
clt 


01 


D -2 


Our notation is A a bk a k b = Akk and akab&k b = °" 2 - Now, we define an entropy density for any 
generic higher curvature theory as: 


S = ^ f (1 + p) Vh d D 2 x, (6) 

J c 

where p(t) contains contribution from the higher curvature terms including ambiguities. For 
general relativity, p = 0. In the stationary limit p will coincide with the Wald expression ([Tj). 
Now from this expression the change in entropy per unit area gives us the following expression 

3 Interested readers are referred to [23 EH for discussions on the Raychaudhuri equation. 
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of generalized expansion 0 , 


e ~f + <4(1 + p) ' 


(7) 


The evolution of 0 is governed by the following equation, 

^_ K0= _ l+p)-a|(l+p) + 0*^ 

TT . d 2 P o d P 

+ tikk H—rry — pttkk — , 

at 2 at 

where we have used the equation of motion Q and inserted eq. (J5]) . We can now write eq, 
in a convenient form 


( 8 ) 


K0 8T[Tkk T Ekk , 
at 


where 


dp 


Ekk — Hkk + ^77 _ pRkk + k a k b 'V a 'V b p 
dt 


* + °k ) (1 + P) ■ 


D — 2 


(9) 

( 10 ) 


We have used 


and 


d kl = k a k b V a V b p + n d ^ 
dt z dt 


d 


Next, consider a situation when a stationary black hole is perturbed by some matter flux obeying 
NEC. The perturbation can be parametrized by some dimensionless parameter e. Note that 
Tkk is linear ((9(e)) in perturbation and so as 9k, a k and dp/dt. Now, to establish a linearized 
second law we ignore higher order terms in eq. ([9]). Then, it is easy to see then eq. (10) reduces 
to, 

Ekk — VfcV kP — pRkk + Hkk ■ (11) 

We have already mentioned that Tkk is of order e, so if the rest of the terms in ([9]) are also 
collectively of higher order, i ,e.,Ekk ~ E>(e 2 ), then we obtain 


— - K0 = - 87 r T kk . (12) 

The above equation implies dO/dt — k 0 < 0, on every slice of the horizon. We assume that 
in the asymptotic future, the horizon again settles down to a stationary state, we must have 
0 —>• 0 in the future. This will imply that 0 must be positive on every slice prior to the future 
and as a result the entropy given by ([ 6 ]) obeys a local increase law. 

Therefore, to establish the linearized second law, we only need to show that the linear order 
terms in E kk exactly cancel each other. Interestingly, this alone will be enough to obtain the 





values of both the coefficients p and q introduced in the entropy functional (16). In the next 
section, we will demonstrate this explicitly for the curvature squared gravity theories. 


4 Linearized second law and fixing JKM ambiguities 

We start with the most general second order higher curvature theory of gravity in five dimen¬ 
sions. The action of such a theory can be expressed as, 

S= lL J {R-2A+ aR 2 + f3 R ab R ab + 7 C gB ) (13) 

where Cgs = R 2 — 4 R 2 ab + R 2 abcd is the Gauss- Bonnet (GB) combination and A is the cosmo¬ 
logical constant. It is reasonable to assume that such a theory admits a stationary black hole 
solution as in the case of GR. We also expect to have a non stationary black hole solution with 
in-falling matter by perturbing this solution. Such a spacetime will be the counterpart of the 
Vaidya solution in general relativity for spherically symmetric case and can be expressed as, 


ds 2 = —f(r, v)dv 2 + 2 dvdr + r 2 dTi\ (14) 

£3 can be any three dimensional space with positive, negative or zero curvature. We want to 
use this solution to investigate the issue of second law of black hole mechanics. Note that, 
the location of the event horizon r — r{v) for this solution can be obtained by by solving the 
following equation, 

, _ dr(u) _ f(r,v) 
dv 2 

with appropriate boundary condition. Note that, (') and (') denote respectively derivative with 
respect to v and r. The null generator of the event horizon is given by k a = {1, f(r, u)/2, 0, 0, 0} 
and the corresponding auxiliary null vector l a = { — 1, 0, 0, 0, 0}. The event horizon has nonzero 
expansion due to the perturbation caused by in falling matter. Next, we will write the entropy 
associated with the horizon as, 


5 


—2?r 



dL 

9 Rabcd 


tabbed P dk@l q&k&l 


(16) 


We consider several choices of the coefficients a, /3 and 7 etc and study the evolution of this 
horizon entropy. The aim is to fix the unknown coefficients p and q for different gravity theories 
by demanding the validity of the linearized second law. 
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4.1 Gauss-Bonnet gravity 


We will first study the local increase law for Gauss-Bonnet (GB) case. This corresponds to the 
choice a = =0 in (13). The action is: 

dPxJ^c 


S = 


16n 


R - 2A + 7 ( RabcdR abcd ~ 4 R ab R ab + R z ) + C 


ab 


(17) 


where we have also introduced a matter sector obeying the NEC. 7 is a coupling constant of 
dimension Mass~ 2 . The equation of motion for this theory is given by, 


With, 


G a b + Rgab + Hab ~ 8nG T a b. 


(18) 


H a b ~ 27 


p p OD P c O P c ^ P _i_ p cde p 

~ Zlt ac lt 5 — Zlt itacbd 1 it a itbcde 


-jQabt-GB- 


(19) 


Next we start with the expression of the entropy functional ([ 6 ]) and evaluate the entropy 
density for the EGB theory. This is given by sum of Wald entropy density plus the ambiguity 
terms, 


S = - J (1 + p) Vhd 3 x, 

(1 + 2 'q(R + AR kl - 2 R km -p9 k 0i - qcr k Vi)) Vhd 3 x , ( 20 ) 

where 2 (R + 4 R ki — 2 R k iki) is the contribution from the Wald construction. Also p = 2 r yp and 
q = 2 7 q. The EGB theory belongs to the general Lovelock class of action functions which give 
rise to quasi linear equation of motions. It has already been shown in CD that black holes in 
all Lovelock theories obey a linearized second law if one uses the JM entropy functional. For 
EGB gravity, the JM entropy functional is the intrinsic Ricci scalar ( 1Z ) of the horizon slice 
and using the null Gauss-Codazzi equation we can cast 1Z as, 



R — R + 4 R k i — 2 R k iki ~ -9k@i + 2cr k ai- (21) 

We will now explicitly show below that in EGB gravity the entropy functional which obeys 


the second law is indeed the JM entropy by fixing the coefficients p and q in (20). 

First, we proceed to evaluate the rhs of eq. 0 to study the second law. Note that, if we use 
the metric (30) with the choice of to be a flat metric then the shear term in ([ 9 ]) will vanish 
identically and q will remain undetermined. This happens because isometries of the metric on 
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a horizon slice essentially coincides with that of a sphere. So we will break the symmetry by 
adding a cross term and the metric on the horizon slice takes the form, 


ds 2 = —f(r, v)dv 2 + 2 dvdr + r 2 (dx 2 + dy 2 + dz 2 + e\h(r, v)dxdy). 


( 22 ) 


We will assume that this shear mode h(r,v) will be balanced by some matter stress tensor still 
obeying the NEC. Also we do not require to find the explicit form of h(r,v) for our analysis. 
We will calculate Ekk order by order in e± and extract the coefficients of the linear order terms 
in e from the evolution equation. Setting those terms to zero will satisfy the linearized second 
law and in the process p and q will be determined. Below we quote expressions for 9's and cr's 
for this solution. 

3 f (r qA 

(23) 


e k = 3 4^d + o(el), 


e,= 


2 r{y) 
3 

r(v) 


+ 0(e 


and 


CW/ — 


e i ( r(v) 2 h'(r , v) 2 — 3 r(v)h(r)h'(r, v ) + 2 h(r, v ) 2 ) f(r, v) 4n 

Mvf + ° (£l) ' 


(24) 


(25) 


Now evaluating E^k and extracting the zeroth order terms in t\ (and linear order in e), we 
get the following equation, 

d 2 f(r,v ) 


d 2 v 


-(4-3 p) = 0. 


(26) 


From 0(e 2 ) terms (which is when h(r,v ) makes its first appearance) we get another equation, 

2(—6 + 6 p + q)h(r, v) 2 + 3(2 — 3 p — q)h(r, v)r(v)h'(r, v) + (2 + q)r(v) 2 h'(r, v) 2 = 0. (27) 

ft is evident that, both of these equations are satisfied if 

4 

P=3 q = - 2. 

This shows that linearized second law fixes all the quadratic ambiguities in entropy functional 


(16) uniquely and it is the JM entropy which obeys the linearized second law for GB black holes! 


Although we have shown this calculation using the 5-dimensional lagrangian for concreteness, 
the same result holds in any dimension D > 4. 


4.2 R? ab theory 

Let us repeat the above analysis for other curvature squared theories. First we take R 2 lb theory 


and we set a = 7 = 0 and (3 > 0 in (13). We will start with the following entropy functional, 
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5 


(28) 


= \J d 3 xVh(l - 2 /3(R k i - p6 k 6i - qa k a t )J. 

When p and q are zero then this reduces to the corresponding Wald entropy for the stationary 
case. Again proceeding as before we will use the linearized second law to fix these coefficients. 
Like the Gauss-Bonnet case we get two independent equations and solving them we get, 



q — 0. 


As argued in ra. the shear part has not contributed to the entropy functional at linear order 
and the other coefficient of the entropy functional exactly matches with the one which has been 
shown to obey the linearized second law in |15j . 


4.3 R 2 theory 


We now consider the case 3 — 7 


0 and a > 0 in (13) and the entropy functional now reads, 


5 = 


1 

4 


d 3 xVh(l + 2 a(R — p9 k 0i — qcr k .ai 


(29) 


It has already been argued in [T3] that f(R) theory obeys a linearized second law and the Wald 
entropy functional itself does the job. It is clear, in this case setting p and q equal to zero will 
correspond to the Wald entropy. If we again repeat the linearized second law analysis we find, 


p = 0 q — 0. 


So this is exactly what is expected from the earlier analysis in 


Therefore, considering all the individual terms in (13) we have demonstrated that for any 


general curvature squared theory we can fix the entropy functional completely using only the 
linearized second law. 


5 Beyond linearized second law 


In this section, we will consider local entropy increase law beyond linear order which means 
we are not allowed to neglect the terms in (J9]) which have been thrown away in the previous 
section. Now we will be considering all the relevant terms in equation (J9]) and will determine 
the criteria such that second law holds non-perturbatively in the coupling. For this, we will 


assume a spherically symmetric Vaidya like solution. This means, the metric in (14) will now 
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have the following form, 

ds 2 = — f{r , v)dv 2 + 2 dvdr + r 2 dVt 2 . (30) 

In fact for EGB gravity one can obtain such a solution just setting the mass in the static 
spherically symmetric Boulware-Deser |2SJ solution to be a function of the advanced time v. 

Beyond linearized level, one has to evaluate the full E k k in the evolution equation (|9]). Since 
the solution is spherically symmetric, shear is identically zero. In fact, for this particular case 
the evolution equation will be of the form, 

^ = -Svr T kk - cel (31) 

Note that we must have dO/dt < 0 to have a local entropy increase law. Also, we have T kk > 0 
by NEC. Now consider a situation where the stationary black hole is perturbed by some matter 
flux and we are examining the second law when the matter has already entered into the black 
hole. In that case, the above evolution equation does not have any contribution from matter 
stress energy tensor and the evolution will be driven solely by the B\ term. So we will have a 
equation of the form, 

f-e = -cc (32) 

In such a situation, if we demand the entropy is increasing, we have to fix the sign of the 
coefficient of term. We evaluate the coefficient in the stationary background and impose 
the condition that overall sign in front of Q\ is negative. This will give us a bound on the 
parameters of the theory under consideration. 


5.1 Gauss-Bonnet Gravity 


We start with the EGB gravity. The event horizon is now a null surface whose equation is 
r = r(v). Calculating the r.h.s. of d9l) for the metric (30) we obtain, 


dQ_ 9/(r,u) 2 

dt ^ 4r(u) 2 


(33) 


where we have identified 

2 = 9/M) 2 
k 4 r{v) 2 

and introduced ( as the coefficient of the d\ terms. The expression of ( is given by, 


c 


1 

3 


l + 2 7 



2 df(r,v) \ 
r(v) dr J 


(34) 
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For 7 = 0 , the coefficient reduces to 1/3 which matches with GR. 1Z is the ricci scalar 
evaluated on the horizon. To satisfy the entropy increase law we now set, 


C > 0 . (35) 

As discussed earlier, we will evaluate ( for different stationary backgrounds and determine 


bounds on the coefficient 7 imposing the condition (35). 


5.1.1 Asymptotically flat case 

Now we consider the Boulware-Deser (BD) [28J black hole as the background, for which, 


/ M 




(36) 


Also, 


r l + 27 = M 


(37) 


determines the location of the horizon. Existence of an event horizon demands r\ > 0. I 11 this 
case horizon topology is a sphere. Now, evaluating f for the above background at the horizon 
r = rft, we get the following inequality 


M 2 + 4My + 36y 2 
4y 2 - M 2 


(38) 


which leads to 


M > 2 17 1 if M > 0 (39) 

M < — 2 |y| if M < 0. (40) 


To understand this better, note that we require M > 2y to avoid the naked singularity of the 
black hole solution for 7 > 0. Thus in this case for a spherically symmetric black hole ( will be 
positive and hence second law will be automatically satisfied. The condition of the validity of 
the second law is same as that for having a regular event horizon. 

Also, for 7 > 0, it is possible to make r h as small as possible by tuning the mass M. But 
when 7 is negative (a situation that appears to be disfavoured by string theory, see [25j, [ 2 Sj and 
references therein), cannot be made arbitrarily small and it would suggest that these black 
holes cannot be formed continuously from a zero temperature set up. Notice that we could 
have reached the conclusion without the second law if M is considered to be positive-however, 
our current argument does not need to make this assumption. Due to this pathology, it would 
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appear that the negative GB coupling case would be ruled out in a theory with no cosmological 
constant. 


5.1.2 AdS case 


Next we consider the 5-dimensional AdS black hole solution for EGB gravity as the background. 
The function /(r) now becomes [17] . 


/m = * + J(i-^-£M?)‘)- <«> 

Here l is the length scale with the cosmological constant (2A = —12 /l 2 ) and k can take values 
0,1 and —1 corresponding to planar, spherical or hyperbolic horizons respectively. The intrinsic 
Ricci scalar on the horizon is 6 k/r\ where, the horizon is at 



(42) 


First we will consider black brane solution for which k — 0 and the horizon 
intrinsic curvature of the codimension two slices of the event horizon vanishes, 
coefficient ( reads, 



is planar. The 
In this case the 

(43) 


Introducing a rescaled coupling A gbI 2 = 2y [THJ [30j we get, 


C — g(l ~ 8A gb ). 


(44) 


Again demanding positivity of ( we get, 

A GB <g, 


(45) 


so to satisfy the second law using S at 0(e 2 ) order we need to impose a bound on the GB 
coupling. Incidentally the same has to be imposed to avoid instabilities in the sound channel 
analysis of the quasinormal mode for dual plasma. It was shown in [30] that when Agb > | 
the Schroedinger potential develops a well which can support unstable quasinormal modes in 
the sound channel. It is interesting to see that the second law knows about this instability. 

Next we consider k — — 1 case. This will give us a black hole with a hyperbolic horizon. We 
will also set ro = 0, which corresponds to the zero mass solution. In this limit we find , 


C — 5\/l — 4 A g _b — 4. 


(46) 
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Demanding ( > 0 we get, 


(47) 


A GB < 


9 

100 ' 


If we take the limit r$/l —» oo in the expression of ( with k = ±1 (i.e. the case when the 
horizon sections becomes planar), we recover the bound on GB coupling A gb for the k = 0 
case. However, this is a weaker bound on A gb and the strongest bound on A gb comes from the 
hyperbolic black hole in the massless limit. This is illustrated in figures la and lb. We have 
checked that for the k — 1 case other thermodynamic considerations (e.g., positive entropy, 
black hole being the correct phase) do not lead to stronger bounds. 



Figure 1: (Colour online) (a) Plot of ( for different ro for k = —1. (b) Plot of ( for different 
r 0 for k = 1. The black dot denotes A gb = 9/100. We have set the 1 = 1. Note that both the 
figures shows that ( will be positive in all the cases provided A gb < 9/100. 


Quite curiously we have recovered the tensor channel causality constraint [20] [18] H9j/] The 
bound on A gb hr this channel leads to the lower bound on the ratio of shear viscosity ( r /) to 
entropy density (s), Anri/s > 16/25 in GB holography [THJ . The vector and scalar channels 
also lead to bounds on the coupling from causality constraints. The strongest bound in GB 
arises from the tensor and scalar channels [19]. Our finding coincides with the tensor channel 
constraint. It will be interesting to see if the other channels can also be reproduced using the 
second law analysis. 

We end this section with one last comment. Instead of putting a bound on the Gauss- 
Bonnet coupling, one can try to modify the entropy functional S itself so that all the 0(e 2 ) 
terms will be canceled. One possible modification is to add a term like Q\Q 2 which vanishes on 
the stationary horizon and boost invariant in the null direction. So this type of higher order 

4 We can repeat the same analysis for higher dimensions using the normalizations in ||19j . Curiously, we 
find the same bound A gb < which is stronger than the causality bounds for d > 5 in [T5j. For flat case 
(k = 0), the bound on GB coupling in d dimensions can be expressed as A gb < 1/(2 (d— 1)). This means for 
any dimension the strongest bound arises when the black holes have hyperbolic horizon (fc = —1) with ro = 0. 
This would also suggest that in the large d limit, ??/s —> 1/47T contrary to 1/87T obtained in [T9 j . 
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ambiguous term can occur in the Wald derivation. We modify S as, 


Smodify = 7 / d 3 xVh( 1 + 2y U + xOltf)- 


( in (43) will be modified as, 




1 - 


47 df(r,v) + 54y f df(r,v ) 


r(v) dr r{vY“ 


dr 


( 48 ) 


(49) 


So we can adjust x such that the ( vanishes and second law automatically holds up to 0(e 2 ). 
One may also can add a term like at this order where V; is the covariant derivative 

with respect to the surface index. But these additional terms will generate further higher 
order terms starting from 0(e 3 )- So one will need to add more terms to cancel them. It might 
possible to do this recursively. However, the coefficients of these terms will also depend on the 
background data. So it will be prudent not to modify the entropy functional. We have put a 
bound on the coupling instead of modifying the entropy functional at 0(e 2 ) order such that 
second law continues to hold up to this order. It will be interesting to investigate what happens 
to this bound on A GB when one goes beyond 0(e 2 ) order. 


5.2 R 2 ab theory 

5.2.1 Asymptotically flat space 


First consider the asymptotically flat black hole with spherically symmetric horizon in 5d. In 
this case, 


C = )( 1 + 7V 

3 V r[v ) 2 




Pf'(r,v)~ 


25 


iPf{r,v) 


(50) 


2r(u) 2 

Now, to extract a bound, we need an explicit form of a spherically symmetric static black hole 
solution in Ricci 2 theory. We are not aware of any such solution except the 5d Schwarzschild 
solution of GR which is also a solution of Ricci 2 theory. So, as an example we take such a 
solution as the background and write, 


/(f) 



(51) 


Then the validity of the second law gives us the condition, 


■r 2 
' h 

34 


>(3. 


(52) 
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When (3 < 0 this is automatically satisfied whenever a horizon exist whereas for (3 > 0, a lower 
bound on the horizon size ensures the validity of the second law. 


Repeating this same analysis for 4d with the metric function f(r,v ) = 1 — (r^/r) we get , 


r 2 
' h 
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(53) 


Therefore, in both cases, there is a minimum horizon radius for the second law to hold when 


(3 > 0. Again as discussed in 5.1.1 minimum ry scenario would appear pathological and may 
be taken as a reason to disfavour the /3 > 0 sign. 


5.2.2 AdS case 

We repeat the same analysis of the previous section for R 2 ab theory but with the background as 
an asymptotically AdS black brane in 5d. The metric function is then given by 

f( r ) = ~pfoo (l - (' r h /r ) 4 ) (54) 

where foo = (yyy:) 2 and it satisfies the equation, 

1-/oo + ^A 2 /^ = 0. (55) 

In this case we obtain, 


c= (l + ^/>.«) + 


Pf'(r,v) 

2 r(v) 


25 Pf(r,v) \ 

6 r(v) 2 ) ’ 


(56) 


Dehne A 2 Z 2 = 2/3 and after evaluating on the horizon with the above metric function as back¬ 
ground, we get 


C = \/9 - 24A 2 - ? 


(57) 


From this, the bound becomes, 


A 2 < 


17 

216' 


(58) 


In this case if we repeat out analysis for a zero mass hyperbolic black hole we obtain a lower 
bound on A 2 . Combining them we get, 


29 x 17 

-A 2 -. 

243 216 


(59) 
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If we repeat the same analysis for 4d and obtain the following bound, A 2 < pp No analog of 
such a bound in the context of AdS/CFT is known in this case. This analysis also suggests 
that the second law bound arising from hyperbolic horizons is not necessarily connected to the 
causality constraints. The reason is that for these theories, in the analysis of [20], we only get 
07 1 > 0. 


5.3 Critical gravity 

We will now consider bounds arising from the second law in quadratic curvature theories with 
a negative cosmological constant, which only involve terms like R a bR ab and R 2 . In particular, 
we will focus on the case of critical gravity [3T1 ,32]. For any theory with only Ricci 2 or R 2 , the 
positive energy condition of Hofman-Maldacena ra does not yield any constraint (see [33] for 
general expressions) except for the positivity of the two point function of stress tensors which 
follows from unitarity. In all the examples below, the positivity of black hole entropy leads to 
the same condition arising from demanding ct > 0 where Ct is the coefficient in the two point 
function of CFT stress tensors. The second law constraint will lead to further conditions as we 
will see. 


5.3.1 D=4 


The action for the 4d version of critical gravity is given by, 

s = -h / d'x^g (r + p + al 2 ^R ab R° b - -R‘ \ ) 


(60) 


In general this theory posses a massive spin-2 mode in addition to the usual massless spin-2 
degree of freedom. It has been shown in [31], at a = —1/2, which is called the critical point, 
the additional massive spin-2 mode becomes massless. Also, at the critical point the entropy 
of the Schwarzschild-AdS black hole becomes zero. 

We will now use the second law to put a bound on the coupling a. Proceeding as before at 
the second order we get, 



2 al 2 f'(r,v) 2 2 2 2 A 

-- + o ' a 1 f ( r w) * 0 a 1 r ( v )f A v) . 

Sr(v) 3 3 / 


(61) 


In this case, we take the background as, 


/M = 




(62) 
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From this we get the bound on the coupling as, 


a < —. 
“ 12 


( 63 ) 


Also, the black hole entropy has to be positive-if it was zero then it must remain zero at all 
times. From this we have the following condition, 


1 + 2a > 0. 


(64) 


So, 


Combining (63) and (65) we get, 


0> --. 


1 i 

- < a < —. 

2 ~ 12 


(65) 


( 66 ) 


Now from the AdS/CFT perspective it is natural to assume the positivity of the coefficient ct 
of two point correlator of boundary stress tensor. That will give exactly the same condition 


as in (65) which is weaker than (66). Repeating the same analysis for a zero mass hyperbolic 


black hole we did not any further bounds. 


5.3.2 D=5 

Next we consider the critical gravity theory in 5d. The action is given by, 

1 


S = 


167r 


J d s x^9 [R+j +c, l 2 (j R ab R“ b - If 


(67) 


The critical point now corresponds to a = —b/27 [31j. Proceeding as before at the second order 
we get, 

(l al 2 f'(r,v ) 5 2 5 ,2 / w ///, A 

C=( v 3 + 2r ^ +-^(*lf(r,v)-—alr(v)f (68) 

In this case, the background solution is taken as, 


nr^ r M-m 


(69) 


where the quantity satisfies the equation, 1 — foo ~ i a f7o — 0- The second law thus leads 
to, 


1 109 

- < a < -. 

3 “ “ 2809 


(70) 
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Again demanding the entropy to be positive we get, 1 + |/oo a > 0, from which a > — A. Thus 
we have, 

5 109 

-< a < -. 

27 “ 2809 

If we repeat our analysis for a zero mass hyperbolic black hole, the resulting bound we obtain 
is weaker than this. 



5.3.3 D=3 

We end this section by exploring the case of NMG theory in 3d. The action is given by, 


S = 


167T 


J d?x\/—g ^R + + ct l 2 ^R a bR a 




The critical point is now at a = —3. For this case, 

7 al 2 f'(r,v) 9 


^ V 1 2 r(v) 

The background solution is taken as, 


+ ^l 2 f"(r,v) - -al 2 r(v)f"(r,v)y 


'M-W 1 -ft) 1 )- 

where the quantity /qq satisfies 1 — foe + a /<T — 0. Then we get the following bound, 


9 

a < — . 
“ 25 


(72) 


(73) 


(74) 


(75) 


Demanding the positivity of the black hole entropy we get, 1 + ^|-^ > 0. This gives, —3 < a < 1. 
Thus we have, 

(76) 


9 

3< “ £ 25 - 


In this case, even if we consider a zero mass hyperbolic black hole we obtain the same bound. 


6 Holographic c-functions from the entropy evolution 
equation 

In this section we will construct a holographic c-function using the FPS entropy functional for 
general curvature squared theories. Here we will again use the evolution equation for 0 ([9]) and 
in the process we will have a geometric derivation for the c-function. 

For unitary Lorentz invariant theory one can define a function in the space of couplings 
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which decreases monotonically from UV to IR, parametrizing the RG flow. It coincides with 
the central charges of the corresponding CFTs at the UV and IR fixed point. In 1 + 1 and 3 + 1 
dimensions this has been established rigorously pRE, 35j. Now from the holographic point view 
this c function plays a pivotal role in understanding nature of RG flow. A proposal for the c 
function for arbitrary theories of gravity has been given in (2T, 37. J3Sj for arbitrary dimensions. 
In this section we will use the Raychaudhuri equation as discussed in the previous sections 
and derive a c function for a general curvature squared theory. For general higher curvature 
theories, we do not expect a c-theorem as there are bound to be problems with unitarity in 
such theories. Thus some condition on the couplings, which presumably only explore a subset 
of conditions leading to unitarty as in [21], is expected. This is what we will find as well. At 
the UV and IR fixed points this function will coincide with the A-type anomaly coefficient for 
the curvature squared theory and it will also have to be a monotonically decreasing quantity 
along the RG flow. To prove the monotonicity we will use Raychaudhuri equation. The basic 
setup is similar as before. 

6.1 Setup 

To derive a c function purely in terms of geometric quantities we will first start with a specific 
example. We will consider domain wall type of geometry in 5 dimensions, the metric for which 
is given below, 

ds 2 = dr 2 + e 2A ^ r \—dT 2 + dx 2 + dy 2 + dz 2 ). (77) 

The spacetime can be foliated by surfaces of codimension-2 at every constant time slice. The 
induced metric on such surfaces is given by, 

dslurface = e 2A(r) ( dx 2 + dy 2 + dz 2 ) . (78) 

Next we will consider a hypersurface orthogonal null congruence. The tangent vector for this 
congruences is given by, 


k a = ( d x ) a = {-e~ A{r \ e~ 2A{ - r \ 0, 0,0}. (79) 

It satisfies k a \7 a k h = 0 and k a k a = 0. We can construct another auxiliary null vector l a as, 

A(r) 

r = (1 2~’2’ 0 ’ 0 ’ 0} (80) 

such that, k a l a = —1. Note that in this case k a and l a are not any kind of horizon null generators 
as compared to the previous sections. They simply correspond to a null congruence that 
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converges along the light sheet projected out of the codimension-2 surface under consideration. 
Also in this section we will use affine parametrization. 


6.2 Holographic c-functions in four derivative gravity 

Next we will start with a guess for the c function for curvature squared theory. The action for 
this theory is given by, 


S = 


167T 


<f'X 


R + an + (3R ab R ab + 7 ( R 1 - 4 R ab R ab + R abcd R abcd ) 


(81) 


The corresponding entropy functional is the FPS functional. 

Sfps = j f d 3 xVh( 1 + p) = j [ d 3 xVh( 1 + 2 aR — j3(R k i - Xd k 0i) +7^ j. 


(82) 


1Z is the 3 dimensional Ricci scalar intrinsic to the surface. Now we will start with the following 
candidate, 


c(A) = 


0 


y/hOi 


(83) 


where as usual, 0 = 0^(1 + p) + ^ and A is the affine parameter. This appears to be a natural 
guess since for the Einstein theory it goes over to the c-function in [23j. Then we have to show 
two things, 


1. c(A) is monotonically decreasing under the flow along the null geodesic congruences. That 
is we have to check the sign of 

2. Further it has to coincide with the correct central charge (A-type anomaly in four dimen¬ 
sions) at the fixed points. 


Now, 


dc( A) 
dX 


1 dG 
1 dG 

TKeldx 


0 d 

(Vhetydx 


(Vhe 


9 


Vhei 


(84) 


using = \/hd k . Next we will replace ^ 
pi). We obtain, 

dc( A) 
dX = 


by (8) with k 


0 and 


24:/r GT kk 

Vhd 4 k 


+ Ekk j 


^ by Raychaudhuri equation 
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where, 


Ekk 


- ((^ 13 >(r) + 15.4'(r) 3 + 16 J 4'(r) J 4"(r))(15 J 4'(r) 4 + 2(^“>(r)+ 

8X(r) 2 ) + 82 ll 3 >(r)A'(r) + 6,4'(r) 2 ,4"(r))) + WA+|±M2/ - 2i(A^(r) 

+ 8 A"(rf) + 60A'(r) 4 (4(10a + 2/3 + 3^)A"{r) - 3) + A'(r) 2 A"(r) 


(256(10a + 2/3 + 3rf)A"(r) - 207) + 8A (3) (r)A'(r)(4(10a + 2/3 + 3^)A"(r) - 15)) 

4e- A ( r) (10a + 2/3 + 37 ) 24 " (r) 

27A'(r) 2 ‘ 


( 86 ) 


Also, 


dc( A) 
dA 


_ 2A< 4 '(r)/T(r) + 8A+r).4"(r) - 16A< 3 >(r).4'(r) 2 
15A'(r) 3 A"(r) + 4$A\r)A"(rf) - M , (f)4 l > 


(87) 


and 


c(A) = 


2(16a + 5 /9)A®(r) + A'(r) ((16a + 5/3 ) (16A"(r) + 15A'(r) 2 ) - 6) 

162A'(r) 4 ' 


( 88 ) 


Notice that in (87) the term proportional to 16a + 5/3 has no possibility of having a definite 
sign and it does not cancel with terms in E^ so it is natural to set the coefficient to zero. Now 
Vhei is positive and T^k using the null energy condition is positive. With 16a + 5/3 = 0, we 
ffiid ^ 

(89) 


d / 4(10a + 2/3 + 3y) \ 

i\\ 27A'(r) / 


Euu — — 
dX 

where we have used A — k a W a . So we will have now, 


dc( A) 2A'KGTkk d / 4(10a + 2(3 + 3y) 

~dX~ Vh6* + dAV 27A'(r) 

Then we can define a effective quantity, 


such that, 


c(A) 


c(A) + 


4(10a + 2/3 + 37 ) 
27A'(r) 


dc{ A) _ 2AnGTkk 
dA “ y/hd 4 k ~ ’ 


(90) 


(91) 


(92) 
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Next evaluating c(A) and multiplying both side by a factor of — ^ we get, 


d /l-4(10a + 2/3 + 37)A'(r) 
dA V A'{rf 

Now in terms of the radial coordinate r this becomes, 


< 0 . 


(. k r d r + k T d r 


1 - 4(10a + 2/3 + &y)A'(r) 
A'{r) 3 


< 0 . 


( 93 ) 


(94) 


From this, 


-2 A(r) A ( l ~ 4 ( 10a + + 3 r r) A 

dr V A'(r) 3 


> 0 . 


(95) 


So , 


1 - 4(10a + 2/3 + 3y )A'(r) 2 
A'(r) 3 


is a monotonically decreasing quantity along the RG flow, where the RG scale is determined 
by the radial coordinate r. At the fixed point the geometry is AdS and hence 


T 

A{r) ee —. 

^ads 

The c(r) coincides (upto an inconsequential overall factor) with the A-type anomaly coefficient 
at UV and 1R (see eg. [33j for expressions in general higher curvature theories). In this 
construction, the starting point of the entropy functional played an important role in reaching 
the condition 16a+ 5/3 = 0. In [2T1I36] , this condition emerged by demanding the existence of a 
“simple” c-function. Our analysis gives another way of looking at the same condition. One may 
wonder if it is possible to come up with a scenario where we did not need to absorb E kk into 
the c-function like what we have done. For example, instead of the HEE functional, perhaps 
one could try to construct a functional that cancels the E kk on the rhs of the Raychaudhuri 
equation. This does not appear to be the case-starting with the Wald entropy functional 
instead of the HEE functional does not naturally lead to the 16a + 5/3 = 0 condition. 


7 Discussion 

In this paper, we considered four derivative gravity and subjected it to the condition that 
the second law of black hole thermodynamics is satisfied at linear and quadratic orders in 
perturbation. We showed that the entropy functional is uniquely fixed (upto second order in 
extrinsic curvatures) and coincides with the holographic entropy functionals. We also derived 
interesting bounds on the coupling which arises at next to leading order in perturbations. We 
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conclude here with a brief discussion on possible future directions. 

• One should consider more generic perturbations to see what happens at linear and sec¬ 
ond order. This is important since the possibility exists that for generic perturbations 
the second law holds only for very specific cases like the Lovelock theories. In [37] , it 
was argued that causality constraints would require adding an infinite set of higher spin 
massive modes in a GB theory for consistency. It will be interesting to see if the second 
law knows about this in the perturbative sense used in the current work. 

• By considering asymptotically flat black holes and the Gauss-Bonnet theory, we found 
that for the negative GB coupling we get a minimum horizon radius for the second law to 
hold. Taking this to be a pathology, this appears to disfavour this sign of the coupling. 
We further saw that the second law appears to know about the sound channel instablity 
in the context of Gauss-Bonnet holography. It will be interesting to see if A > —1/8 which 
corresponds to the scalar channel instability [T8] can also be seen from the second law. 

• One should consider how our analysis will change if there is matter coupling to the higher 
curvature terms. In this case, it is not clear how to define the null energy condition. It 
seems natural that the matter couplings will get constrained by demanding the second 
law. 

• Our analysis can be extended to general theories such as the quasi-topological theories 
[38] , In order to get a sufficient number of conditions to fix all the extrinsic curvature 
terms, one will need to turn on generic perturbations. 

• In case of gravitational Chern-Simons terms, it will be interesting to see if and how 
the validity of second law can fix the entropy functionals [39]. Since the rhs of the 
Raychaudhuri equation does not know about topological terms, the linearized analysis 
will not be able to capture the effect of such terms. One may need to consider topology 
changing processes like black hole mergers [25] to see the effect of such terms. 

• It will be interesting to compare our analysis with what arises from similar considerations 
in the fluid-gravity correspondence [TO]. Demanding a positive entropy current for higher 
curvature theories will lead to constraints which are presumably going to be similar to 
what we have found. 

• The connection between the second law constraints and the positivity of relative entropy 
[IT] in the context of holography can also be explored. We found that when we consider 
the second law for topological black holes with zero mass, we recover the bound from the 
tensor channel causality constraint. Why is this happening? A plausible explanation is 
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as follows. The topological black holes with zero mass are just AdS spaces written in 
different co-ordinates and the finite Wald entropy of these black holes can be interpreted 
as an entanglement entropy across a sphere 1211 (2! in the dual field theory. Thus the 
second law in this context may be related to the positivity of relative entropy due to a 
time dependent perturbation corresponding to the infalling matter we have considered. 
Relative entropy is expected to be positive for a unitary field theory. Hence there appears 
to be an interesting interplay between bulk causality and field theory unitarity. It will be 
interesting to make this connection more precise. It will also be interesting to understand 
the bounds in the context of c-theorems using the second law for causal horizons as 
advocated recently in [42]. 
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